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Summary 

The mechanism by which the noise generated at the blade passing frequency by a propeller is 
altered when the propeller axis is at an angle of attack to the freestream is examined herein. The 
measured noise field is distinctly non axially symmetric under such conditions with far field sound 
pressure levels both diminished and increased relative to the axially symmetric values produced with 
the propeller at zero angle of attack. Attempts have been made to explain this non axially symmetric 
sound field based on the unsteady ("once per rev") loading experienced by the propeller blades when 
the propeller axis is at non zero angle of attack. A calculation based on this notion appears to greatly 
underestimate the measured azimuthal asymmetry of noise for high tip speed, highly loaded propell- 
ers. A new mechanism is proposed herein; namely, that at angle of attack, there is a non axially sym- 
metric modulation of the radiative efficiency of the steady loading and thickness noise which is the pri- 
mary cause of the non axiaily symmetric sound field at angle of attack for high tip speed, heavily 
loaded propellers with a large number of blades. A calculation of this effect to first order in the 
crossflow Mach number (component of freestream Mach number normal to the propeller axis) is car- 
ried out and shows much better agreement with measured noise data on the angle of attack effect. 
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Introduction 


In recent years, there has been a renewed interest in high tip speed turboprop systems for com- 
mercial aircraft applications. Many of the practical applications envisage counter rotating propeller 
(CRP) driven aircraft though careful experimental evaluations (both aerodynamic and acoustic) have 
been made of both single rotation propellers (SRP) and CRP’s. 

In the present study, we address a particular aspect of the noise radiated at the blade passing fre- 
quency (BPF) by a propeller. In an ideally uniform stream with the propeller axis coincident with the 
direction of the free stream, aerodynamic sound at the BPF produced by the propeller blades is due 
to the steady air loading (drag and lift) carried by the blades and also due to the thickness of the 
blades. A characteristic of this sound field is that it is axially symmetric in magnitude with the axis of 
symmetry being the propeller axis. The first theoretical analysis of steady loading noise was given by 
Gutin (1936) with Garrick & Watkins (1954) contributing the extension of Gutin’s work to the case of 
a propeller in forward flight. Much progress has been made since these studies in terms of our ability 
to predict steady loading/thickness related noise. The two principal avenues of improvement in 
recent years have been the use of computational fluid mechanics techniques to predict the steady 
pressure distribution on the blade surfaces and the recognition in the acoustic theory of source non- 
compactness effects. 

When the propeller is installed in an aircraft, clearly it does not operate in the environment of a 
uniform stream whose direction is coincident with the propeller axis and the sound field at the BPF 
can be quite different from the previously described ideal case. These effects are described as "instal- 
lation effects". An elementary example of the installation effect is that under takeoff/approach condi- 
tions, the direction of the propeller axis deviates from the direction of forward flight leading to the 
"angle of attack" effect. Recent experimental studies of the angle of attack effect on propeller BPF 
noise include those of Block (1984) and Woodward (1987a) for SRP’s and of Woodward (1987b) for 
CRP’s. Paradoxically, as will be discussed below, this most elementary installation effect has proved 
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quite difficult to understand. 

The characteristics of the measured non axially symmetric sound field at angle of attack can be 
described in terms of the polar coordinate system shown in figure 1. The BPF sound field (which is 
axially symmetric at zero angle of attack) roughly speaking is found to be diminished in the upper 
quadrants on the right hand half of Figure 1 (i.e. for 0< <f> < v) and increased in the lower quadrant (tt 
< <f> < 2 rr) with the minima and maxima occurring approximately at <j> = jt/ 2 and <f> = 3 tt/ 2 relative to 
the zero angle of attack case. For highly loaded, high tip speed propellers the changes can be quite 
large. In the vicinity of the plane of rotation the changes can be as large as 1 dB per degree of angle 
of attack. All prior attempts to explain this noise change due to angle of attack have proceeded as fol- 
lows. The BPF noise of the propeller due to steady loading and thickness is calculated using the usual 
methods for a propeller in uniform axial flow. This sound field due to steady loading/thickness which 
is axisymmetric in amplitude but not in phase is assumed to be unchanged with the propeller operating 
at angle of attack. Due to angle of attack, clearly unsteady loading at "one per revolution" will be 
experienced by the propeller blades which will cause additional noise at the BPF. When this new 
sound field is added to that due to the steady loading with proper consideration being given to the 
phasing between the two sources, a resultant non axisymmetric sound field is obtained. 

As will be shown in the next section, for lightly loaded propellers especially with few blades and 
spinning at low tip speeds, this procedure does seem to yield results in agreement with data. For 
highly loaded, high tip speed propellers especially with larger number of blades, this procedure yields 
gross underestimates of the measured departure from the axisymmetrical situation due to angle of 
attack. A new mechanism is pointed out which negates the notion that the sound field due to steady 
loading/thickness is unchanged by angle of attack. A first order (in angle of attack) estimation of this 
mechanism is carried out which appears to greatly improve the theory-data comparisons. Indeed the 
suggestion that emerges from this analysis is that for highly loaded, high tip speed propellers with a 
large number of blades, the angle of attack effect on BPF propeller noise is largely unrelated to 
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unsteady loading. 


The notation used in this report is indicated in Appendix 2. 
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Steady Loading and Unsteady Loading Noise 

In the work described herein, the acoustic source distributions are assumed to be acoustically 
compact in the chordwise direction. In other words, a line source model is employed. Also in this 
section the only mean flow effect accounted for is a uniform axial flow corresponding to the flight 
velocity. At time "t = 0," let the location of this line source be as shown in Figure 4, and let the sweep 
be as shown in Figure 5. Thus, (r)" and "Ax (r)” define the radial lean and axial sweep of the 
blades. The following results were derived for the purposes of this work by the use of axial Fourier 
transforms and the method of stationary phase, but they are in agreement with results given by Hawk- 
ings and Lowson (1974) and Hanson (1984), and hence, only the final result of the analysis is given. 

The contribution to p Q " from the blade element lying between "r" and "r + dr" is given by: 


dp 0 _ 

dr 


jnB 2 n j 

2nRc (l-M F cosff) ±nB 


nBflr sin 0 
c ( \-M f cos0) 


exp [ j [ ±nB W-f) + nBa i( r )- nBQt 

1 J (I-Mfcosff) 2 J 

' fzo cosfl - cf^_ + j nBp 0 cQhl 1 (1) 

(I-Mpcosff) + Hr (i-MpCosfff 

In Equation (1), the far field noise is expressed in terms of R, 6 (the coordinates of the observer 

relative to the "retarded" propeller disk location). The relationship between these coordinates and 

R , 0 (the coordinates of the observer relative to the current propeller disk location) is shown in Fig- 
0 c 

ures 2 and 3. In Equation ( 1), wherever two signs appear, the upper sign refers to counterclockwise 
rotor rotation, and the lower sign to clockwise rotation (forward looking aft). 

Now consider the noise due to unsteady loading. Let the unsteady loading per blade per unit 
spanwise length in a frame of reference fixed with a rotating blade be expressed in the form. 
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(2a) 


and 


Unsteady z component 


E f A cos(kQt) 

k =1 


Unsteady <f> component = E cos(£Qf) (2b) 

Both f zk and f 4> k are P ositlve functions of "r". Also, let /? zk (r) and ^(r) denote the azimuthal 
angles in the coordinate system of Figure 1 where these unsteady components achieve their maximum 
positive values. Then, unsteady loading noise from blade elements lying between "r" and "r + dr" con- 
tributes to the far field pressure a quantity dp k /dr, given by the sum of: 


d Pk) 


dr 


and 


df>k_ 

dr 


where 


\dp k ] 

jnB 2 0 

nBOr sin 9 

dr 

< > 

4nRc (UMfCOse) ' (±nB+lc) 

c ( 1-A/pCos 9) 


exp 


) 


&(./?+ Ax cos 9) 
(1-MfCos 9) 


+ ( ±nB+k ) 



+ nBcti (r) - nBDt 


exp (-jkfa) facosS c exp (-jkfa) / # (nB ±k) 
(l-A/ f cos 9) nBfir 


( 3 ) 


and 


' d P k 

jnB 2 0 

nBOr sin 9 

dr 

< > 

" 4 xRc (I-Mpcos 6) 

c(l-Mpcos 9) 


exp 


/ 


k(R +Ar cos 9) 
( 1-MfCos 9) 


+ ( ±nB-k ) 



+ nBa x (r) -nBftt 


exp (jkfa) f^cos 9 - c exp (jkfa) fa 
(1-M/rCos 9) nB Or 


(nB + k) 


( 4 ) 
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The total far field acoustic pressure can now be computed as: 


P = / 


E 

dr * = 1 


(dPk) 


dr 


00 
+ E 




dr 


dr 


( 5 ) 


where r. and r denote the inner and outer radii of the propeller disk, respectively, and where p 
h t 

is a complex number whose amplitude gives the amplitude of the far field pressure in the nth har- 
monic of blade passing frequency. 


The above formulation could, perhaps, be expressed more concisely with more extensive use of 
complex notation. However, it was felt preferable to restrict consideration to positive n and k and, 
also, to consider both clockwise and counterclockwise rotating propellers. These interests, and the 
desire to be accurate concerning phase relations, have led to the above formulation. 


The theory needed to compute the steady axial and tangential forces for a SRP is based on that of 
Glauert (1963). The formulation is based on incompressible flow theory and the inputs needed are: 
relative radius (r/R) of interest, solidity (chord/transverse spacing) at this radius of interest, number 
of blades, specification of the section lift and drag coefficients as functions of angle-of-attack (includ- 
ing the orientation of the reference line with reference to which the angle-of-attack is specified), and 
the speed ratio, A = U/Q R The version used for the work described herein incorporates the Prandtl 

tip loss factor. 

For the CRP case, an extension of Glauert’s theory by Lock (1941) with inclusion of an approxi- 
mate accounting of slipstream contraction is used. 

In order to ensure that lift and drag coefficients as functions of angle of attack are realistic, these 
relations were deduced from performance maps for the relevant propellers by inverting the 
Glauert/Lock theories at a representative span location such as 70% span. This however entails 
several additional assumptions such as that for the flow per unit area and work addition per unit mass 
of fluid are uniform. We assume that the coefficient relations deduced in this manner for the 70% 
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span location are valid for all spanwise locations. 

To determine the unsteady forces, in this study, we have used a quasi-steady approximation based 
on the Glauert/Lock theories to deduce the unsteady forces. The angle of attack induces a low fre- 
quency modulation of the blade loading in a frame of reference fixed with the rotor blades and hence 
a quasi-steady approximation should be adequate. The quasi-steady approximation for SRP’s and 
CRP’s is outlined in Appendix 1. 

A fairly complete theory data comparison with data from Block (1984) based on the theory of this 
section can be found in Whitfield (1989). In Block (1984), the tests were carried out on a four bladed 
SRP designated as an SR2. Some sample results from Whitfield (1989) are reproduced herein. 

Figure 6 shows results for the axisymmetric sound field at zero angle of attack at a tip speed of 
443 fps. In Figure 7 a theory-data comparison is shown for the change in noise due to angle of attack 
in the plane of rotation. Figures 8 and 9 show similar results for the same propeller at a higher tip 
speed of 745 fps. It is a characteristic of the theory outlined in this section that, in the plane of rota- 
tion, the acoustic factor for BPF noise determining the ratio of unsteady loading to steady loading 
effects scales as "Jb Since this ratio diminishes as "M t " (wheel tip Mach number) 

increases through subsonic values, it is not surprising that smaller angle of attack effects are predicted 
for a higher tip speed (Figure 9 versus Figure 7). 

The results in Figures 7 and 9 are encouraging considering the data scatter though it is arguable 
that a systematic underprediction of the magnitude of the noise decrease is evident in Figure 9. 

In Figures 10(a) and 10(b), we show the results of applying the theory of the present section to 
angle of attack effects measured by Woodward (1987b). These measurements are on a 11x9 bladed 
CRP spinning at approximately 800 fps in a Mach 0.2 airstream. The theory-data comparisons for the 
angle of attack effect on the BPF noise of the front rotor indicate a severe underprediction of the 
measured effect. 
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The key difference between the propellers tested by Block (1984) and Woodward (1987b) is that 
the latter propellers operate with substantially higher steady loading (lift coefficients are higher by a 
factor of 2-3) and the number of blades involved in the Woodward experiments (11 and 9) is much 
greater than the four blades used in the Block experiment. 

In the course of analyzing the disappointing results of Figure 10(a) and 10(b) (an exhaustive 
theory-data comparison based on the results of the present section was carried out with all the BPF 
data available in the Woodward experiments with results typically as shown in Figure 10) the author 
was made aware by colleagues at the NASA Lewis Research Center that presently unpublished 
research has yielded similar gross underprediction of the measured azimuthal variation of BPF noise 
with angle of attack. Estimates of the unsteady force levels (derived herein by a quasi-steady approxi- 
mation) are of course of interest in the context of forced vibration of the propeller blades and conver- 
sations with colleagues in that discipline indicate that the unsteady force estimates provided by the 
present procedure for the angle of attack effect are reasonable. It also seems necessary to contend 
with the observation that the theory of the present section appears to perform adequately in explain- 
ing the Block (1984) experiments with lightly loaded, four bladed propellers. 

The experience cited above has led to a re-examination of a basic premise of the present develop- 
ment; namely, that the steady loading/ thickness noise itself is unaffected by angle of attack. As the 
discussion of the next section will show, this premise appears incorrect. An approximate accounting 
of the effect of angle of attack on steady loading/thickness noise is carried and a very substantial 
improvement is now evident in terms of the Woodward (1987b) studies. The development of the next 
section will also elucidate why this effect is relatively unimportant for lightly loaded, low tip speed pro- 
pellers with a few number of blades as were involved in the Block (1984) experiments. 
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The physical mechanism which causes the steady loading/thickness noise of a propeller to be fun- 
damentally non axially symmetric (i.e. in both amplitude and phase) when the propeller is at angle of 
attack is quite simple and can be explained in terms of Figure 11. 


In this figure we show a steady source executing a circular motion as would be characteristic of steady 
sources rotating with the blade. Also shown in the figure is the uniform upwash provided by the cross 
flow due to angle of attack. Now the radiative efficiency of a source depends on the relative velocity 
between the source and the fluid surrounding it and as Figure 11 shows this relative velocity (and 
hence the radiative efficiency) is modulated in a non axisymmetric fashion due to the cross flow. This 
effect has been neglected in prior analyses. 

To capture this effect, we proceed as follows. Including the cross flow (designated as "My" in Fig- 
ure 12), the acoustic pressure due to steady loading satisfies: 

Pxx Cl -Ml) +Pyy (\-M 2 y ) + p n + 2jkMj> x + 2jkM yPy - 2M x M yPxy + k 2 p = (— + — --) 

dx t d<f> 

where k = uj/c, all quantities have a time dependence as exp (-jwt), "c” is the speed of sound, and 
My are axial and cross flow Mach numbers, "J^, f^" are the force terms and j = v^T and p xx denotes 
cPp/dx 2 , etc. In the following, M x will be denoted by M. 

Let 


fx = fxo exp 6 (*) 6 ( r-a ) 

f* = f*o exp ( jm<f > ) 6(x) S (r-a) 


where <f> is the azimuthal coordinate. In the present study, it has been possible to solve the above 
problem only to O(My). 


Define the axial Fourier transform pair: 
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oo 


P = f p exp (-jox)dx 

-OO 

and 


p = — / P exp ( jox ) da 
1 jL 


Then to O (My), P satisfies: 


Pyy+Pn+ 2 IjMyPy (k-aM) + [ (k-aM) 2 - a 2 ] P = j of„ + (™ )/* 
Now let P = P exp jjyA/ y (fc-oA/| 

Then to O (My), P satisfies: 

Pyy + Pzi + [ {k-aM) 2 - a 2 ] P = ; [a/ TO + ra/^/a] 


5 (r-a)ei‘ 


im<f> 


8{r-a)ei m + 


1 + M y (k-aM) a ( - ~ - ) 


Let 


7 = V" (k-aM) 2 -o? 

Then 

P = y [(<**) /*> + ™/*>] <*P (fc-aiWj 
[**"♦./* (to) (r) + M ^ k ^ )a * >(, ” +1 » 

/*♦! (to) flj&i (T) - e j{m-Wj m l (t a)/f£> x ( 7 r)] 

We can now set x = -R cos 5 , r = R sin 0 and y = R sin 0 sin ^ and carry out an asymptotic 
c c c c c c 
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evaluation of "p" for large R £ . The full details of this asymptotic evaluation are not necessary and it 
suffices to note the following. Firstly in the evaluation of the "m" azimuthal mode in the above there 


are O(My) shifts in the point of stationary phase (relative to the case My=0). This shift of the point 
of stationary phase is illustrated in Figure 13, in terms of the resulting new relation between current 
and retarded coordinates. It is reiterated that Figure 13 is valid only to O(M^). Note that the relation 
between current and retarded coordinates itself now is non axially symmetric i.e. it depends on 
However this detail is far less significant than the appearance of new azimuthal modes of order (m-1) 
and (m+ 1) in the equation above. For the "m-1" and "m+ 1" modes we can use the point of stationary 
phase as for "My=0" since these are O(My) contributions and hence relative to the "m" mode, these 
contribute as: 


m-1 


- exp [-M ~ V 2 1 M y {ka) j 
2(1 -Mcosff) m ' 1 


ka si n 6 

(1 -M cos 6) 


m + 1 


exp \j(^/2)]M y (ka) 

2(1 -Mcosff) m+1 


■v 

ka sin 6 
(1 -M cos 9) 


where the relation between "0 " and retarded angular coordinate "0" is as shown in Figure 3. Clearly 
the "m-1" term contributes much more significantly than the "m +1" term. A virtually identical treat- 
ment applies to thickness noise. Note that "ka = nBM t " where "n" is harmonic of BPF noise, "B = 
number of blades" and = wheel tip Mach number" at radius "a". 

The modifications to equation (1) to O(My) will now be written out. In order to clearly motivate 
this development, this procedure will be carried out step-by-step. Firstly we reconsider and rewrite 
equation (1). Introduce the following symbols. 
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Symbol 

Value 

Significance 

prmlt 

jnB 2 il/(2nC) 

premultiplier 

camp 

1/(1-M cos#) 

convective 
amplification factor 

arg 

rtBf2rsin(#)(camp)/c 

argument of 
Bessel function 

pht 


phase 


+ nBai(r)-nBCtt]] 

term 

tht 

jr \Bp 0 cDhl 

thickness 

term 


With the above symbols, the present form of equation (1) can be written as: 

-f- = < P’ mU } <plu) 

dr R 

+ J ±nB (arg) + (tht)(campf J ±nB (arg)) 

Ur 

When effects on steady loading/thickness noise to O(M^) are calculated, the following new symbols 
are introduced. 

i. aat (for "angle of attack" term): 

( aat ) = (nB QrM y {camp )/2c) 

(exp(j(^rt/2))J ±nB + l ( arg ) - 
exp (-](<(* r/2)) J ±nB -\{arg)) 

ii. = My sin <f> (component of cross flow Mach number in ”<f>” direction) 

iii. A# = M y<j , (M-cos #) /( 1-M cos 6) (change in emission angle to O(M^) due to cross flow for fixed 
R c- 9 c> 

iv. cos (ip) = cos# - A# sin 6 
sin(^) = sin# + A# cos # 

(V is the new emission angle to 0(M ) for fixed R £ # c ) 
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V. R n = R(1-M^/(1-A/COS0)) 


(R n is the new distance of observer from center of retarded propeller to O(My) 
*c> 

vi. campn (for new "convective amplification factor") = 1/(1 - A/cosi/> + A/^sinfl) 

vii. argn (for new argument of Bessel functions) = nBOr(sinrp+ M y<t ,)(campn ) /c 

Then the result analogous to equation (1) including the crossflow effect to O(M^) i 

dPo 


dr 


= (prmlt ) (pht ) \f zo (cos xjj{campn) 2 


j ±nB { ar gn)/R n + (, aat ) cos (6) ( camp) 2 jR ) 


- c fjo {campn )f gn) (aat)(camp) . 

Hr ( R n R } 

+ (tht){(ccimpnf J ±nB (argn)/R„ + (aat)(campf /R)\ 


for fixed R 
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Revised Theory-Data Comparisons and Discussion 


The results presented in this section use the revised acoustic theory for steady loading/thickness 
noise and the old acoustic theory (wherein the only mean flow effect dealt with is a uniform, axial 
flow) for the unsteady loading/thickness noise. This approach is consistent with the idea that at 
present only an 0(M y ) acoustic theory has been developed where "My" is the crossflow Mach number 
due to angle of attack. All "ff' coordinates referred to in Figures 14-21 are current coordinates (as 
opposed to retarded coordinates). All predictions are for the fundamental BPF tone of the relevant 
propeller and are for the noise change due to angle of attack. 

In Figures 14-15, using the new acoustic theory, we show theory-data comparisons for the Block 
(1984) data. It can be seen that the comparisons are at least as good as shown in Figures 7 and 9. 
The new effect discussed in the previous section scales as "nBM t " and pertains only to steady 
loading/thickness noise. Thus it is not surprising that the new effect is not very consequential for a 
li gh tly loaded, four bladed propeller such as the SR-2 used in the Block (1984) experiments. 

In Figures 16-21 we show the theory data comparisons for the CRP data of Woodward (1987b). 
Several comments concerning these theory-data comparisons are noted below. 

1. In cases where the measured T is not close to a multiple of 10 degrees, predictions for multi- 
ples of 10 degrees bracketing the measured T are both compared with the data. 

2. Figures 16(a), (b) may be compared to Figures 10(a), (b) (the data shown in these figures are 
the same). A dramatic improvement is now evident with the revised acoustic theory in terms of 
the substantial under prediction of the old acoustic theory of the angle of attack effect. 

3. Figures 17 and 19 indicate that angle of attack effects on the rear rotor noise are still 
significantly underpredicted. The rear rotor in Woodward (1987) has fewer blades (9 versus 
11) than the front rotor and hence the predicted magnitude of the new effect (recall the "nBM t 
scaling law for the new effect) is correspondingly smaller. Also as far as unsteady effects are 
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concerned, the steady and quasi-steady CRP theory developed herein predicts a "smoothing" 
effect of the front rotor in terms of the effect of inflow distortions on the rear rotor and hence 
the unsteady effects on the rear rotor are also predicted to be smaller than on the front rotor. 
At any rate the reasons for underprediction of effects on the rear rotor are not clear. 

4. Finally if we examine figures 18(a), (b), 19 (as a typical case), it would appear that the maxima 
and minima of the data appear at larger ”<f>' with respect to the theory for the front rotor and 
vice-versa for the rear rotor. For the CRP used in Woodward (1987b), the front rotor spins in 
the direction of decreasing (clockwise forward looking aft). It may be surmised that when an 
elongated centerbody such as used in Woodward (1987b) is placed at an angle of attack to the 
freestream and the propellers are placed aft of the maximum diameter location of the center- 
body, the propeller blades could experience a once per rev distortion of incoming axial velocity 
such that this velocity is diminished on top and increased on the bottom. In Figures 20(a), (b) 
and 21, a calculation was carried out assuming the presence of such an axial velocity distortion 
of magnitude 10% of the freestream velocity. Comparing these results with figures 18(a), (b) 
and 19, clearly a shift is obtained of the theory which is in the right direction relative to the data. 
The purpose of this exercise is to indicate that small distortions of axial velocity additional to the 
tangential velocity distortion created by angle of attack can cause significant shifts of the predic- 
tions in the "0" direction even if the magnitudes are not much affected. 

A complete version of theory-data comparisons with all of the data obtained by Woodward 
( 1987b) will be published in a NASA contractor’s report to follow. 
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Concluding remarks 


In this study an approach based on steady and quasi-steady aerodynamics of a propeller has been 
developed to predict both the steady loading/ thickness BPF noise and the unsteady loading/thickness 
BPF noise developed by a propeller operating at angle of attack. Originally a conventional acoustic 
theory was used in which the only mean flow effect accounted for is a uniform axial flow. Application 
of this approach to data of Woodward ( 1987b) indicated a severe underestimation of the magnitude of 
the effect of angle of attack on BPF noise from a CRP. A new idea has been developed in this paper 
based on the notion that the small crossflow due to angle of attack (designated in Mach number as 
"My") actually profoundly filters the steady loading/thickness BPF noise (causing it to be azimuthally 
non uniform). This effect may be far more important than unsteady loading/thickness noise for 
heavily loaded, high tip Mach number propellers with a large number of blades. This effect has been 
calculated only to O(My) in this study. Revised theory-data comparisons indicate that the underesti- 
mation problem is substantially resolved especially for BPF noise from the front rotor though 
underestimation of the magnitude of noise decreases due to angle of attack for BPF noise from the 
rear rotor is still in evidence. Since the O(My) effect has itself proved to be large, it is suggested that 
research is needed to calculate the acoustic impact of the crossflow more fully. 
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Appendix 1: Quasi Steady Aerodynamics of Propellers 


Single Rotation Propellers 

Let u(r,<£), v(r,<£) denote externally contributed nonuniform axial and tangential velocities incident 
on the propeller, "u" is taken positive in the downstream direction and "v" is taken positive in the 
direction of propeller rotation. 


At any radius, the quasi-steady forces (thrust and tangential force) are obtained in the functional 
form: 


A) 

where "ft" denotes the propeller rotational velocity in radians per second, and "A" the speed ratio 
"U/ftR," where U is the advance velocity of the propeller. 

The change in the force, F, due to u and v can be obtained as follows. 


In the quasi-steady approximation, changes in ft and A due to u and v, will be: 


5ft = - 


6 A = A 


+ — 
U ftr 


The resulting change in F can be written: 


6F = 2ft/(A)5ft + ft 2 



6 A 


-2F + A 

dF 

— + 

f 

dF 


dX 

ftr 

d\ 

K 


u 

nR 


Thus, if u, v are nondimensionalized by introducing: 
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u ' = and v 
HR 


v 

OR 


we have: 


6F = [-2F + A — v ' + ^ « 
3A r d\ 


Now, dF/dX can be approximated by a central difference formula: 


dF ~ F(A+e)-F(A-e) 
dX ~ 2e 


where V is a small increment. 
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Counter Rotation Propelle rs 


Now, in general, the extension of the steady single rotation propeller theory to the counter rotation 
case will yield for quantities (Q p ) associated with the front rotor a functional relation of type: 

Qf = C F n}fr (A ,S) 

where "Cp" is a suitable constant of proportionality at a given radius. Likewise, for quantities (Qp) 
associated with the rear rotor we may write: 

Qr = c R nlf R (a ,s) 

We now wish to compute the variation 6Q F , 6Q R , due to the extra (u,v) velocity field. We note the fol- 
lowing intermediate relations: 

SC1 F = -v/r, Sn R = v/r 
SX = A [H/t/+v/(rfV) | 

SS = -S [v/(rfV)+v/(ril/{)| 

= -5(l+S)v/(KV) 

-•S.(l+S). A. i-'-)(v/U) 


Thus: 
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Qf p 

2n F sn F f F + n % (— ■ ) sx + n 2 F (-^-) 


= CfQp 


A A ( + 

A ' 1 dX ) ' K U ) 


i2 

( as'-] 


&s| 


aA'~ ' 'as 

-2. A. (f Ifr H- A 2 . (f ) . (^) -5 . (1 + S) . A. (^) (^r) 
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Likewise: 
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2 . Pr . 5fl/? • /}{ + ^ Or ( ^ ) • ^1 
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i? A 
( as /w J 
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These are the relations that are used to determine the fluctuation quantities in the quasi-steady 
approximation for the counter rotation case. Note that partial derivatives with respect to both "A" and 
S (advance ratio and speed ratio) are now required. These are evaluated by examining steady solu- 
tions for the following five pairs of points in the (A,S) space: 

(A,S), (A+0.01,5), (A-0.01,5), (A,S+0.01), and (A,S-0.01). 
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Appendix 2: Notation 


B Number of propeller blades 

c Ambient speed of sound 

Cj^: Constants of proportionality associated with front and rear rotors at a given radius, r 

f zo Steady (thrust) force per unit spanwise length exerted by each blade in the z direction (Fig- 

ure 2); f ZQ is a function of radius V 

ffr Steady (tangential) force per unit spanwise length exerted by each blade in the "</>" direction 

(Figure 1); note that is a function of radius and also that, in view of Figure 1, if the 
sense of rotor rotation is clockwise (forward looking aft), will be negative 

h Representative airfoil thickness as a function of radius 

j V-1 

J Bessel function of the first kind 

k Wave number, nBCl/c 

1 Airfoil chord as a function of radius 

Mp Flight Mach number 

n Harmonic of blade passing frequency of interest (n > 1) 

p Q Far field acoustic pressure due to steady loading and thickness 

r Radius of interest 

R Distance of observer from center of retarded propeller disk location (reference is made to 

Figure 3) 
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Distance of observer from center of current propeller disk location (Figure 3) 


R Propeller radius 

P 

S Rotor speed ratio, i.e., S =ttf/Q R 

t Time 

U: Forward flight speed 

u: Perturbation axial velocity component (positive downstream) 

v: Perturbation tangential velocity (positive corresponds to "v" in same sense as direction of 

rotation of forward rotor) 

4, Azimuthal angle as defined in Figure 1 (note that <f> = 0, 180 represent the horizontal direc- 

tion) 

A: Advance ratio based on front rotor, i.e., A =U/ (RCl F ) 

Ambient air density 

9 "Emission" angle for propeller noise 

$ c "Observer" angle for propeller noise 

n Angular velocity of propeller in radians per second 

n F ,n R : Speeds of rotat ion of forward and rear rotor 
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Figure 4. - Definition sketch for radial lean. 
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gure 5. - Definition sketch for axial sweep. 
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16(a), except ©predicted “ 110 
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o 

^r 
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Same as figure 16(a), except rear rotor BPF and 6 measure d = ^predicted 
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Figure 18(a). I I x I) b laded 2 It. diameter CRP; 75% radius pitch angles 35 
axial flow Mach number 0.2; 0 degree ang 1 e of-a t tack ; symbols are measured 
Predicted based on revised theory; ©measured 88 H)4°, ©predicted = 100°; fn 
Woodward (1087b); increase/decrease of noise due to ang le of-at tack . 
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redicted ~ 




47 


Figure 1‘). - Same as figure 18(a), except rear rotor BPF and ©measured = ©predicted - 90 
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Figure 2U(a). - Same as figure 18(a), (front rotor BPF , ©measured ~ 104 > ^predicted - 
angle-of-at lack) , except that a 10% once per rev axial velocity distortion is assumed 
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Figure 21. - Same as figure 20(a), except for rear rotor BPF and ©measured = ^predicted 
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